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We report on the successful experimental generation of electron bunches with ramped current 
profiles. The technique relies on impressing nonlinear correlations in the longitudinal phase space 
using a superconducing radiofrequency linear accelerator operating at two frequencies and a current- 
enhancing dispersive section. The produced ~ 700-MeV bunches have peak currents of the order of a 
kilo- Ampere. Data taken for various accelerator settings demonstrate the versatility of the method 
and in particular its ability to produce current profiles that have a quasi-linear dependency on 
the longitudinal (temporal) coordinate. The measured bunch parameters are shown, via numerical 
simulations, to produce gigavolt-per-meter peak accelerating electric fields with transformer ratios 
larger than 2 in dielectric-lined waveguides. 

PACS numbers: 29.27.-a, 41.85.-p, 41.75.Fr 



Electron acceleration is a rapidly-advancing field of sci- 
entific research with widespread applications in industry 
and medicine. Producing and accelerating high-quality 
electron bunches within very compact footprints is a chal- 
lenging task that will most probably use advanced accel- 
eration methods. These techniques can be categorized 
into laser-driven [1-3] and charged-particle-beam-driven 
methods [4-7]. In the latter scheme, a popular config- 
uration consists of a "drive" electron bunch with suit- 
able parameters propagating through a high-impedance 
structure or plasma medium thereby inducing an electro- 
magnetic wake. A following "witness" electron bunch, 
properly delayed, can be accelerated by these wakefields. 

CoUinear beam-driven acceleration techniques have 
demonstrated accelerating fields in excess of GV/m [9, 
10]. The fundamental wakefield theorem [S] limits the 
transformer ratio - the maximum accelerating wakefield 
over the decelerating field experienced by the driving 
bunch - to 2 for bunches with symmetric current pro- 
files. Tailored bunches with asymmetric , e.g. a linearly- 
ramped, current profiles can lead to transformer ratio 
> 2 [11]. To date, there has been a small number of 
techniques capable of producing linearly-ramped elec- 
tron bunches. A successful experiment demonstrated the 
production of 50-A ramped electron bunches using sex- 
tupole magnets located in a dispersive section [12] to im- 
part nonlinear correlation in the longitudinal phase space 
(LPS). Unfortunately, the method introduces coupling 
between the longitudinal and transverse degrees of free- 
dom which ultimately affects the transverse brightness of 
the drive and witness bunches. 

In this Letter we present an alternative technique that 
uses a radiofrequency (rf) linear accelerator (linac) oper- 
ating at two frequencies. It has long been recognized that 
linacs operating at multiple frequencies could be used to 
correct for LPS distortions and improve the final peak 



current [13, 14]. We show analytically and demonstrate 
experimentally how a two frequency linac could be op- 
erated to tailor the nonlinear correlations in the LPS 
thereby providing control over the current profile. 

We first elaborate the proposed method using a ID- 
IV single-particle model of the LPS dynamics and take 
an electron with coordinates (^,(5) where z refers to the 
longitudinal position of the electron with respect to the 
bunch barycenter (in our convention z > corresponds 
to the head of the bunch) and 6 = p/{p) — 1 is the frac- 
tional momentum spread {p is the electron's momentum 
and (p) the average momentum of the bunch). Consid- 
ering a photo-emission electron source, the LPS coor- 
dinates downstream are {zo,6o = uoZq + b^z^ + 0{zq)) 
where oq, and 6o are constants that depend on the 
bunch charge and operating parameters of the electron 
source. For sake of simplicity we limit our model to sec- 
ond order in zq and 5q. Next, we examine the accel- 
eration through a linac operating at the frequencies /i 
and /„ = n/i with total accelerating voltage V{z) = 
Vi cos{kiz + ipi) + Vn cos(fc„z -f tpn) whcrc Vi^„ and ipi,n 
are respectively the accelerating voltages and operating 
phases of the two linac sections, and fci^„ = 27r/i^„/c. 
In our convention, when the phases between the linac 
sections and the electron bunch are = the bunch 
energy gain is maximum (this is refer to as on-crest op- 
eration). Under the assumption ki^nZo <C 1 and neglect- 
ing non-relativistic effects, the electron's LPS coordinate 
downstream of the linac are (z/ = zq,Si = a/zg + ^/^o) 
whcrc a/ = ao - e(A:iVi sin (/Si + fc„K sin v3„)/i?/, bi = 
bo — e{kfVi cosiySi + k'^Vn cos(p„)/(2i?/) with e being the 
electronic charge and Ei the beam's average energy down- 
stream of the linac. Finally, we study the passage of 
the bunch through an achromatic current-enhancing dis- 
persive section [henceforth referred to as "bunch com- 
pressor" (BC)]. The LPS dynamics through a BC is ap- 
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proximated by the transformation Zf = R^^Si + T^ee^f 
where i?5g (also referred to as longitudinal dispersion), 
and Tsgg are the coefficients of the Taylor expansion of 
the transfer map {zf\6i) of the BC. Therefore the final 
position is given as function of the initial coordinates 
following Zf = QfZo + bfZQ with a/ = 1 + aiR^Q and 
bf = biR^e + afT^QQ. Taking the initial current to follow 
the Gaussian distribution /o(zo) = IoGxp[—Zq/{2(tIq)] 
(where Iq is the initial peak current), and invoking the 
charge conservation If{zf)dz = lQ{zQ)dzQ gives the final 
current distribution If{zf ) = /o/A-^/^(2;j) exp[— (a/ + 
A^^'^{zj))^/{8bj<jlQ)]e[A{zf)] where A{zf) = a}+AbfZf 
and Q{) is the Heavisidc function. The latter current 
distribution does not include the effect of the initial un- 
correlated fractional momentum spread ct^q. The final 
current, taking into account ag^, is given by the convolu- 
tion Ifizf) = fdz~flf{z~f)cxp[-{zf - i>)V(2a2)] where 
(Ttj = R^eO'gQ. The final current shape is controlled via 
ttf and bf and can be tailored to follow a linear ramp as 
demonstrated in Fig. 1. 
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FIG. 1: (color online) Analytically-computed current profiles 
for several values of b/ for fixed a/ = 2.5 (a) and for several 
values of a/ with bf — 0.7 (b). The numbers in (a) [resp. (b)] 
are the values of &/ [resp. a/]; for all the cases (t„ — 0.05. 

The experiment described in this Letter was performed 
at the Free-electron LASer in Hamburg (FLASH) facil- 
ity [15]. In the FLASH accelerator, diagrammed in Fig. 2, 
the electron bunches are generated via photoemission 
from a cesium tclluride photocathode located on the back 
plate of a 1-1-1/2 cell normal-conducting rf cavity oper- 
ating at 1.3 GHz on the TMqio 7r-mode (rf gun). The 
bunch is then accelerated in a 1.3-GHz and 3.9-GHz su- 
perconducting accelerating modules (respectively ACCl 
and ACC39) before passing through a bunch compres- 
sor (BCl). The ACC39 3rd-harmonic module was in- 
stalled to nominally correct for nonlinear distortions in 
the LPS and enhance the final peak current of the elec- 
tron bunch [!()]. Downstream of BCl, the bunch is ac- 
celerated and can be further compressed in BC2. A last 
acceleration stage (ACC4/5/6/7) brings the beam to its 
final energy (maximum of ~ 1.2 GeV). The beam's direc- 
tion is then horizontally translated using a dispersionless 
section referred to as dogleg beamline (DLB). Nominally, 
the beam is sent to a string of undulators to produce ul- 
traviolet light via the self- amplified stimulated emission 
free-electron laser (FEL) process. For our experiment, 



the bunches were instead vertically sheared by a 2.856- 
GHz transverse deflecting structure (TDS) operating on 
the TMiiQ-like mode and horizontally bent by a down- 
stream spectrometer [17]. Consequently the transverse 
density measured on the downstream Cerium-doped Yt- 
trium Aluminum Garnet (Cc:YAG) scintillating screen is 
representative of the LPS density distribution. The hori- 
zontal and vertical coordinates at the Ce:YAG screen are 
respectively Xs — rjSp, where ry ~ 0.75 m is the horizontal 
dispersion function, and ys — nzp where k ~ 20 is the 
vertical shearing factor and {zftSf) refers to the LPS 
coordinate upstream of the TDS. The exact values of 77 
and K are experimentally determined via a beam-based 
calibration procedure. 



TABLE I: Settings of accelerator subsystems relevant to the 
LPS dynamics used in the experiment and simulations. 



parameter 


symbol 


value 


unit 


ACCl voltage 




[140-157] 


MV 


ACCl phase 




[-10,10] 


deg 


ACC39 voltage 


V3 


[13,21] 


MV 


ACC39 phase 




[160-180] 


deg 


ACC2/3 voltage 


Vl,2-3 


311 


MV 


ACC2/3 phase 


Vl,2-3 





deg 


ACC4/5/6/7 voltage 


K,4-7 


233.9 


MV 


ACC4/5/6/7 phase 







deg 


BCl longitudinal dispersion 




~ 170 


mm 


BC2 longitudinal dispersion 




~ 15 


mm 


Single-bunch charge 


Q 


0.5 


nC 


Bunch energy 


E 


~ 690 


MeV 



The accelerator parameters settings are gathered in 
Tab. I. The nominal settings of BC2 were altered to 
reduce its longitudinal dispersion i?^;^ and the ACC2/3 
and ACC4/5/6/7 accelerating modules were operated on 
crest. Such settings insure that the BC2 and the DBL 
sections do not significantly affect the LPS beam dynam- 
ics. Therefore the measured current profile is represen- 
tative of the profile downstream of BCl. 
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FIG. 2: (color online) Diagram of the FLASH facility. Only 
components affecting the longitudinal phase space beam 
(LPS) dynamics of the bunches are shown. The acronyms 
ACC, BC, and DBL stand respectively for accelerating mod- 
ules, bunch compressors, and dogleg beamline (the blue rect- 
angles represent dipole magnets). The transverse deflecting 
structure (TDS), spectrometer and Ce:YAG screen compose 
the LPS diagnostics. 

In order to validate the simple analytical model 
described above, numerical simulations of the LPS beam 
dynamics were carried using a multi-particle model. 
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The simulations also enable the investigation of possible 
detrimental effects resulting from collective effects such 
as longitudinal space charge (LSC) and beam self inter- 
action via coherent synchrotron radiation (CSR) [1>^]. In 
these simulations, the beam dynamics in the rf-gun was 
modeled with the particle-in-cell program ASTRA [19] 
and the obtained distribution was subsequently tracked 
in the accelerating modules using a ID-IV program 
that incorporates a one-dimensional model of the LSC. 
The program CSRTRACK [20], which self-consistently 
simulates CSR effects, was used to model the beam 
dynamics in the BCl, and BC2 sections. An example 
of simulated LPS distributions and associated current 
profiles computed for different settings of AC CI and 
ACC39 parameters appear in Fig. 3. The results indicate 
that the production ramped bunches is possible despite 
the intricate LPS structures developing due to the 
collective effects and higher-order nonlinear effects not 
included in our analytical model. The simulations also 
confirm that the current profile upstream of the TDS 
(as measured by the LPS diagnostics) is representative 
of the one downstream of BCl. 
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FIG. 3: (color online) Simulated LPS distribution [(a) and 
(b)] with associated current profile downstream of BCl (solid 
blue trace) and DBL (dash red trace) [(c) and (d)] . The set of 
plots [(a), (c)] and [(b), (d)] correspond to different (Vi,3,v'i,3) 
settings. 



Figure 4 displays examples of measured LPS distribu- 
tions with associated current profiles obtained for dif- 
ferent settings of ACCl and ACC39. As predicted, the 
observed current profiles are asymmetric and can be tai- 
lored to be ramped with the head of the bunch (z > 0) 
having less charge than the tail; sec Fig. 4 (b-d). The lat- 
ter feature is in contrast with the nominal compression 
case at FLASH where the LPS distortion usually results 
in a low-charge trailing population as seen in Fig. 4 (a). 

We now quantify the performance of the pro- 
duced current profiles to enhance beam-driven accel- 
eration techniques by considering a drive bunch in- 
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FIG. 4: (color online) Snapshots of the measured longitudi- 
nal phase spaces (left column) and associated current pro- 
files (right column) for different settings of the ACCl and 
ACC39 accelerating modules. The values (Vi, tpi; V3, i^s) [in 
(MV,°,MV,°)] are: (150.5, 6.1; 20.7, 3.8), (156.7, 3.8; 20.8, 
168.2), (155.6, 3.6; 20.6, 166.7), and (156.8, 4.3; 20.7, 167.7) 
for respectively case (a), (b), (c), and (d). 



jccted in a cylindrical-symmetric dielectric-lined waveg- 
uide (DLW) [G] . The DLW consists of a hollow dielectric 
cylinder with inner and outer radii a and b. The cylin- 
der is taken to be diamond (relative electric permittivity 
e,. = 5.7); and its outer surface is contacted with a per- 
fect conductor; see Fig. 5 (a). The measured current 
profiles are numerically convolved with the Green's func- 
tion associated to the monopole mode to yield the axial 
electric field [21]. These semi-analytical calculations were 
benchmarked against finite-difference time-domain elec- 
tromagnetic simulations executed with VORPAL [22]. The 
transformer ratio is numerically inferred asTZ = / E- \ 
where E_ (resp. Ej^) is the decelerating (resp. acceler- 
ating) axial electric field within (resp. behind) the elec- 
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FIG. 5: (color online) Cylindrical-symmetric dielectric-loaded 
waveguide considered (a) and axial wakefield produced by the 
current profile shown in Fig. 4 (c) for (a, b) = (20,60) /im. 




transformer ratio a (urn) 

FIG. 6: (color online) Simulated transformer ratio versus peak 
accelerating field (a) for the four measured current profiles 
(displayed as different colors with label corresponding to cases 
shown in Fig. 4). Transformer ratio (false color map) as a 
function of the DLW inner radius a and dielectric layer thick- 
ness b — a with corresponding \E'^\ shown as isoclines with 
values quoted in MV/m for case (c) of Fig. 4. 



tron bunch; see Fig. 5 (b). The achieved TZ and 
values as the structure geometry is varied are shown in 
Fig. 6. As a e [20, 300] fj.m and b e a + [20, 300] /xm are 
varied the vifavelengths of the excited vifakefield modes 
change. The simulations show that profiles (b) and (c) 
of Fig. 4 can yield values of 7?. > 2. A possible con- 
figuration with (a, b) — (20, 60) /im, results in 7?. ~ 5.8 
with i?+ ~ 0.75 GV/m; see corresponding wake in Fig. 5 
(b). Such high-field with transformer ratio significantly 
higher than 2 and driven by bunches produced in a su- 
perconducting linac could pave the way toward compact 
high-repetition-rate short- wavelength FELs [23]. 

Finally, the proposed technique could be adapted to 
non-ultrarclativistic energies using a two- (or multi-) 
frequency version of the velocity-bunching scheme [24]. 
Such an implementation would circumvent the use of a 
BC and would therefore be immune to GSR effects. 

In summary we proposed and experimentally demon- 
strated a simple method for shaping the current profile 
of relativistic electron bunches. The technique could be 
further refined by, e.g., including several harmonic fre- 
quencies. 
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